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Abstract 



Qh! We present a detailed study of photon production in hadronic events in electron- 

positron annihilation at LEP energies. We show that estimates of the inclusive pho- 
^ . ton spectrum using the quark-to-photon fragmentation function determined using the 

ALEPH 'photon' -|- 1 jet data agree well with the observations of the OPAL collabo- 
ration. This agreement shows that the photon fragmentation function determined in 
^ ■ this way can be used for inclusive observables. We also compare next-to-leading order 

and beyond leading logarithm predictions obtained using the numerically resummed 
solutions of the fragmentation function evolution equation of Bourhis, Fontannaz and 
Guillet and Gliick, Reya and Vogt with the data. Moreover, in order to check the 
general behaviour of the fragmentation function, we consider an analytic series expan- 
sion in the strong coupling. We find that the parameterizations are inaccurate at large 
X values. While the OPAL data is in broad agreement with estimates based on any 
of these approaches, the ALEPH data prefers the resummed BFG parameterization. 
Finally, there is some ambiguity as to whether the fragmentation function is treated as 
0{a) or 0{a/as)- We show that at present this ambiguity affects mainly the prediction 
for the 'photon' -|- 1 jet rate at large z. 



1 



1 Introduction 

The production of hard photons in hadronic processes provides an important testing ground 
for QCD. For example, direct photon production in pp colhsions is used to extract information 
on the gluon content of the proton, while the presence of photons in the final state represents 
an important background source in many searches for new physics. A good understanding 
of direct photon production within the context of the Standard model is therefore essential. 

Photons produced in hadronic collisions can have two possible origins: the d^reci radiation 
of a photon off a primary quark and the fragmentation of a parton i into a photon carrying 
a sizeable fraction x, of the parton energy. While the former direct process constitutes a 
short- distance effect and can be calculated within the framework of the Standard Model, 
the latter is primarily a long distance process. It is described by the process-independent 
parton-to-photon fragmentation function Di^^{x, fip), which cannot be calculated using per- 
turbative methods but which must be determined from experimental data. The evolution 
of Di^.y{x, fip) with the factorization scale fip can however be determined by perturbative 
QCD. 

The most promising environment for a determination of the quark-to-photon fragmenta- 
tion function Dq^^{x, fip) is the study of photon production in electron-positron annihila- 
tion into hadrons. Such measurements have been recently presented by the ALEPH and 
OPAL collaborations at CERN. Both measurements differ not only in the experimental 
observable studied to determine the quark-to-photon fragmentation function, but also in 
the theoretical framework used to match direct and fragmentation contributions onto each 
other. This makes a comparison of both experiments difficult, and has also given rise to the 
speculation that the ALEPH and OPAL analyzes do not probe the same quantity 0, |^. 

The analysis performed by ALEPH is based on the study of two jet events in which 
one of the jets contains a highly energetic photon. Correspondingly, the fraction of the jet 
momentum carried by the photon within the 'photon' jet, z, is greater than 0.7. These events 
were defined by the application of the Durham jet clustering algorithm |Q to both hadronic 
and electromagnetic clusters. A comparison between the observed rate and the 0{a) [|l], |^ 
or 0{aas) theoretical estimates yielded a first determination of the quark-to-photon 
fragmentation function accurate at leading and next-to-leading order. The theoretical basis 
on which the measurement of the 'photon' -|- 1 jet rate relies, is an explicit counting of powers 
of the strong coupling in both the direct and the fragmentation contributions, and where 
no resummation of In fi'jp is performed. We shall refer to this theoretical framework as the 
fixed order approach. 

More recently, the OPAL collaboration has measured the inclusive photon distribution 
for final state photons with energies as small as 10 GeV. This corresponds to the photon 
carrying a fraction of the beam momentum, x-y, to be as low as 0.2. They have compared 
their results with the model-dependent predictions of p, ^ and found a reasonable agreement 
in both cases. These model predictions are based on a resummation of the logarithms of the 
factorization scale hf and naturally associate an inverse power of with all fragmentation 
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contributions. This resummation procedure is the conventional approach. 

Ahhough, the quark-to-photon fragmentation function is a purely non perturbative ob- 
ject, its evolution with the factorization scale can be described within perturbative QCD. 
In the fixed order approach, a byproduct of the perturbative calculation of the photon pro- 
duction cross section is an evolution equation which is accurate at a given order. This has 
been presented in some detail in , and we recall the main features in Section |^. In the 
conventional approach adopted by 0, |^ though, the quark-to-photon fragmentation function 
satisfies the well-known DGLAP p| all order evolution equation and, as in the fixed order 
approach, this solution has a non-perturbative ingredient which can only be measured. We 
discuss the main steps in the formal derivation of this solution in the conventional formalism 
in Section |^. 

A priori, it may indeed not be entirely clear that the fragmentation function extracted 
from a measurement of the 'photon' + 1 jet rate and the fragmentation function arising in 
the inclusive cross section are the same. In fact, the variables used in the definitions of these 
rates, z for the photon rate and for the inclusive rate are generally not equal to each 
other. As explained in Section |], these variables coincide however in the purely 'collinear' 
region where the genuine non-perturbative effects arise so we can clarify this issue and affirm 
that these two fragmentation functions are indeed equal to each other. Information gained 
about the quark-to-photon fragmentation function from the analysis of one observable can 
(and ought to) be used to predict other photon cross sections. 

So far the ALEPH collaboration have compared their data with the theoretical predic- 
tions obtained in a fixed order formalism while the OPAL collaboration have concentrated 
on the theoretical predictions from the conventional framework. However, because the frag- 
mentation functions are universal, it is possible to evaluate the 'photon' -|- 1 jet rate and 
the inclusive rate obtained in either theoretical framework and to compare these predictions 
with either the ALEPH or the OPAL data. To perform such a comparison, and to see if the 
data prefers one approach, is the purpose of this paper. 

More precisely the paper is organized as follows. In Section ^, we discuss how the 
quark-to-photon fragmentation function and the one-photon production cross section are 
determined in a fixed order approach and present the results obtained for the inclusive rate 
in this approach. In Section]^, we describe how the quark-to-photon fragmentation function is 
determined in a conventional approach. In particular we describe how the leading-logarithmic 
(LL) and beyond-leading logarithmic (BLL) solutions are determined and compare their 
numerical parameterizations with analytically expanded expressions of the quark-to-photon 
fragmentation function. Section H contains a detailed presentation of four possible different 
approaches to evaluate the one-photon production cross sections, one of these being closely 
related to the fixed order approach described in Section The possible definitions of the 
non-perturbative input in either the schemes adopted by Gliick, Reya and Vogt (GRV) in 
or by Bourhis, Fontannaz and Guillet (BFG) in 0, are also discussed in this context. 
Moreover we study the behaviour of the parameterizations given in |^, |^ in the large x region. 
In Section |^, we present results obtained following any of these four approaches and using 
either the GRV or BFG schemes for the inclusive and 'photon' + 1 jet cross sections and 
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compare them with the OPAL and ALEPH data. Finally, Section ^ contains our conclusions. 

2 Inclusive photon production in the fixed order ap- 
proach 

Let us first consider the general structure of the single photon production cross section, fully 
differential in all quantities, 

da = da^ + ^ da„ ® Z^f^^. (2.1) 

a 

In this equation which is valid for any single photon production cross section there are two 
contributions. First 'prompt' photon production arises when the photon is produced directly 
in the hard interaction. Second, the longer distance fragmentation process occurs when one 
of the partons produced in the hard interaction fragments into a photon and transfers a 
fraction of the parent parton momentum to the photon. Each type of parton, a, contributes 
according to the process independent parton-to-photon fragmentation functions Da^-y and 
the sum runs over all partons. At the order we are interested in, the gluon fragmentation 
functions will generally be small and can be neglected, as explained below. The sum therefore 
runs over all active quark and antiquark flavours. Moreover, by charge conjugation, quark 
and antiquark fragmentation functions are equal. 

The individual terms in eq. (|2.1| ) may be divergent and are denoted by hatted quantities. 
However, through the introduction of a factorization scale these terms can be reorganized 
and the physical cross section can be written in terms of finite (but factorization scale 
dependent) quantities, 

dcr = dcr^(/ii7) + ^daa® Da^^,{^^F)■ (2.2) 

a 

These two process specific contributions will be defined differently for the 'photon' + 
1 jet and inclusive photon cross sections. Furthermore, as we shall see in Section ^ these will 
be defined differently depending which approach is used to evaluate these photon production 
cross sections. 

We are primarily interested in the expression of the inclusive photon production eval- 
uated at fixed order up to 0{aas)- In this context, the hard cross sections as well as the 
non-perturbative fragmentation functions have to be considered at most up to 0{aas)- Al- 
though the fragmentation functions are non-perturbative, we can nominally assign a power 
of coupling constants, based on counting the couplings necessary to radiate a photon. Since 
the photon couples directly to the quark, Dq^^ is naively of 0{a.) while the gluon can only 
couple to the photon via a quark and Dg_,^ can be considered to be of 0{aas). This simplis- 
tic argument is supported by models of the fragmentation function 0, |^ which suggest that 
gluon fragmentation is a much smaller effect than quark fragmentation. In e^e~ annihila- 
tion, the production of gluons is suppressed by a power of as compared to the production of 
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(a) 



(b) 




Figure 1: Parton level subprocesses contributing to single photon production at 0{aas). 



quarks. As a result, the contribution from gluon fragmentation is of 0{aal) and should be 
neglected in a consistent evaluation of this rate at fixed order up to 0{aas). 

In this fixed order approach the inclusive photon production cross section can be ex- 
pressed as a function of the fraction of the beam energy carried by the photon — 2E^/ -y/s. 
At 0{a) it takes the following form, 



1 dcj^o 



E [[^J Cf{x^,iXF) + D^^^{x^,ix^)]^ , (2.3) 



(7o dx^ 

while at 0{a.a^ it is given by. 



(To dx^ oq dx^ 



+ (£) E ( (^) (^7, y^F) + ® D,_.,{x„ ^^,)^ . (2.4) 



Here, (To is the two-particle cross section, A^^ stands for the number of flavours. We have 
moreover represented the convolution in a compact way, 

J^jA{t)B(^^=A^B{x). (2.5) 



The hard scattering coefficient functions appearing in these equations are defined as 
follows. C^"-* is the coefficient function corresponding to the lowest order process e^e~ qq^. 
It is defined after the leading quark-photon singularity has been subtracted and factorized 
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in the bare quark-to-photon fragmentation function. In the MS scheme it is given by 



In 



III 



(2.6) 



where Pjl^^{x) is the e — part of the the lowest order splitting function in (4 — 2e)- 
dimensions [HI, 



-Pg_^^(x) 



1 



'l-x) 



ex 



X 



(2.7) 



The (finite) next-to-leading order coefficient function C^^-* can be obtained numerically after 
the next-to-leading quark-photon singularity has been subtracted. More precisely, C^^-* is 
obtained after summing all contributions which are independent of Dg^^{x^, fip) arising 
from the Feynman diagrams depicted in Fig. ^ together. A detailed description of the 
evaluation of C^^-* in the case of the 'photon' -|- 1 jet cross section has been given in fTO 



The next-to-leading coefficient function appropriate for the inclusive photon production can 
moreover be straightforwardly obtained from the next-to-leading order coefficient function 
relevant for the 'photon' + 1 jet cross section by integrating the jet-specific variables over 
the complete phase space. 

The coefficient function C^°) is the finite part associated with the sum of real and virtual 
gluon contributions to the process e~^e~ — > qq. It is straightforward to evaluate, and can be 
found for example in H], 



ci'Hx) 



1+x^ 



XI 



-5{l-x)\\n 



+ 



-1-2 In X 



1+x^ 

1 — X 



'ln(l 



x) 



1 — X 



-TT 
J 

3 
' 2 



6(1 



- X) 
3 

X 

2 



,(2.8) 



where Cp the Casimir operator is given by, Cp = {N^ — 1)/2N. This color factor is also 
implicitly included in the next-to-leading order coefficient functions C^^-* defined above. 



As motivated in 0, within the fixed order approach the process- independent quark-to- 
photon fragmentation function Dq^^{x, fip) appearing in eq. (|2.3|) and eq. ( ^^ respectively 
at leading and next-to-leading order, satisfies an exact (up to the order under consideration) 
evolution equation. At next-to-leading order (C(aas)) this equation reads. 



27r 



ae„ 



2n 



^ ] p{i) 
2n 



X] 



+ {^)P^%®Dg^,{x,fiF). (2.9) 
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Pq^q and Pq^y are respectively the lowest order quark-to-quark and the next-to-leading 
order quark-to-photon universal splitting functions 

1+x^ 



11, 12 



Pi%i^) 



Pi%i^) 



Cf 
Cf 



\i-x}+ 

1 9 
1 — X 

2 2 



-5{l~x) 



-8 + -X ) In X + 2x ln(l 



+ ln^(l -x) +41nxln(l -x) +8Li2(l 



X 



X — -TT 



-X In X 
2 



P^,is^) 



(2.10) 



(2.11) 
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The quark-to-plioton fragmentation function satisfying the evolution equation reads, 




s 



in(4)n%®^."^7(^'/^o). (2.12) 



2itJ \fi^ 



This solution has some interesting properties. First, it is exact at the order of the cal- 
culation i.e. 0{aas), and yields no terms of higher orders. The inclusive rate with this 
solution implemented is therefore independent of the choice of the factorization scale hf- 
The exact lowest order (0{a)) evolution equation and its solution D'^^^{x, hf) are natu- 
rally contained in eq. (|2.9| ) and eq. ( p. 12 ) respectively; they can be obtained by dropping all 



terms proportional to in these equations. 

In eq. (|2.12|) , all a priori unknown non-perturbative contributions associated with the 
fragmentation function are contained in Dgl^^{x,fio) which has to be determined from the 
data. This non-perturbative input has been extracted from the ALEPH 'photon' + 1 jet 
data |ll| for ?/cut = 0.06 and z > 0.7. At lowest order, we have |1|, 



D7Ji''\x,fio) = [j^j [-P!,%{x) ln(l - xf - 13.26) , (2.13) 
with /io = 0.14 GeV while at next-to- leading order 0, 

D^^!;'^<'\x,fio) = (^^) (-P(^,(x)ln(l - xf + 20.8 (1 - x) - 11.0?) , (2.14) 

where fiQ = 0.64 GeV and for a,(M^) = 0.124. 

However, in the 'photon' + 1 jet data, the process independent fragmentation function 
is extracted as a function of z, the fraction of the 'photon' jet momentum carried by the 
photon. This is in general different from the variable relevant for the inclusive rate which is 
x^, the fraction of the beam momentum carried by the photon. To see this, let us consider 
the lowest order process e~^e~ qq'j, where the photon is emitted by the quark. For this 
process the two variables x^ and z are defined as follows, 

_ 2-5'7 _ 1 _ _ _ _|_ _ _ -^7 _ Uq-y ~^ Vq-y /r^ i r\ 

where and Eq are respectively the energies carried by the photon and the quark in the 
event and the dimensionless invariants Uij = [pi +PjY/s. Over most of phase space, the 
two variables are clearly different and results derived for the 'photon' + 1 jet cross section 
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Figure 2: The inclusive photon energy distribution normahzed to the hadronic cross section 
as measured by the OPAL collaboration compared with the LO {0{a)) and NLO {0{aas)) 
calculations including the quark-to-photon fragmentation function determined at the appro- 
priate order from the ALEPH 'photon' -|- 1 jet data. 



should in principle not be transferable to the inclusive rate. However, the non perturbative 
fragmentation effects are associated with the emission of a photon collinear to the quark 
and arise in the boundary region of the phase space where i/q^ 0. This is precisely where 
the definition of the two variables z and coincide. Furthermore, the collinear photon and 
quark region of the inclusive three parton phase space, is the same as the colhnear phase 
space restricted to the 'photon' jet region, so that the contribution to the inclusive rate or 
to 'photon' -|- 1 jet rate from the collinear region is the same [0, |], 

--^ = [j^j {P^%i^) In ((1 - x)x) + x) = (2.16) 

where x can stand either for x^ or for z. As a consequence, the quark-to-photon fragmenta- 
tion function obtained in the context of the calculation of the 'photon' -|- 1 jet rate is process 
independent, and can be used directly to estimate the inclusive photon rate. This statement 
is in contrast with claims made in the literature P, ||. 

We therefore show our predictions for the inclusive photon energy distribution at both 
0{a) and 0{aas) in Fig. ^ using the fragmentation function extracted from the 'photon' 
-|- 1 jet data. For comparison we also show the measurements of the OPAL Collaboration 
0. We note that OPAL was able to identify photons with energies as little as 10 GeV, 
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Figure 3: The 0{a) inclusive photon energy spectrum (sohd) broken down into photon + 
1 jet (dashed) and photon + 2 jet (short-dashed) contributions for z > 0.2 and for ?/cut = 0.06. 



corresponding to > 0.2. We see good agreement with the data for both the leading and 
next-to-leading order theoretical predictions. It is also apparent that the next-to-leading 
order corrections to the inclusive rate are of reasonable size indicating that the results ob- 
tained following this fixed order approach are perturbatively stable. The agreement between 
our predictions and the OPAL data is quite remarkable for another reason. The phase space 
relevant for the OPAL data by far exceeds that used to determine the fragmentation function 
from the ALEPH 'photon' + 1 jet data. Let us quantify this statement by examining the 
size of the different contributions entering in the inclusive rate at lowest order. 

At lowest order, at most three jets can be formed with one of them being denoted as 
'photon' jet if it contains a photon. To define 'photon' -|- n jet configurations, we use 
the Durham algorithm |Q. Fig. ^ shows the relative importance for the inclusive rate of 
the 'photon' + 1 jet and 'photon' + 2 jet cross sections for z > 0.2 with ?/cut = 0.06. The 
dominant contribution comes from 'photon' + 1 jet events, precisely those used to determine 
the fragmentation function. However, the ALEPH data used in the fit lies entirely above 
Ey = 32 GeV (corresponding to 2; > 0.7), and the inclusive photon rate for smaller 
should be viewed as a prediction within the fixed order approach. At the very least, the 
good agreement with the inclusive photon data provides some vindication for the functional 
form of the fragmentation function used in fitting the high z ALEPH data. 

In order to quantify the uncertainty of our theoretical prediction for the inclusive rate, 
we have reevaluated it for the two extreme values allowed for /iq as determined by ALEPH 



9 



8*10-" I 

■T • OPAL data 

7*10"'' - ' inclusive 

inciusive - upper fio 




°10 15 20 25 30 35 40 45 

(GeV) 



Figure 4: The 0{a) inclusive photon energy spectrum calculated using the range of uncer- 
tainty on fio obtained from the ALEPH fit to the 'photon' + 1 jet data at 0{a). 



/xo = O.Ut'.il (2.17) 

The results of this calculation together with the inclusive rate for /iq being equal to its central 
value are depicted in Fig. ^. From this figure it appears clearly that the inclusive rate is 
sensitive to the choice of the initial scale /iq but also that the present OPAL data are too 
poor to enable us to reduce the range of uncertainty of /iq. 

The OPAL Collaboration has compared their results with the estimates of [0, p!3|] . 
They find reasonable agreement when the factorization scale /Xi? is chosen to be equal to 
Mz, although the present data were unable to discriminate between the models. In the 
following sections of this paper, we shall describe the different ways with which these results 
were obtained in some detail. As we will see, a common feature of all these model estimates is 
the resummation of leading (a" In""*"^ fip) and subleading (a" In" ^p) logarithms of the mass 
factorization scale fip to all orders in the strong coupling constant. A priori therefore, it 
might seem surprising that the inclusive rate obtained in a fixed order approach as described 
in this section and as depicted in Fig. § and the inclusive rate evaluated in a fundamentally 
different approach where the ln(/i|.) are resummed yield equally good predictions for this 
rate when compared to the OPAL data. In the remainder of this paper, we shall examine 
more closely these two different approaches and one of the major aims of this study will be 
to understand why these two different formalisms yield similar results for the inclusive rate. 
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3 The conventional determination of Dq^^{x, fip) 

In the first part of tliis section we sliall describe liow tlie quark-to-plioton fragmentation 
function is obtained in the so-called conventional approach. In particular we shall describe 
how the leading logarithmic (LL) and beyond-leading logarithmic (BLL) fragmentation func- 
tions are defined within this approach. Analytic expansions of these functions will be given 
in the second part of this section. 



3.1 The conventional approach 

In the conventional approach, the parton-to-photon fragmentation function Di^^ satisfies 
an all order inhomogeneous evolution equation p|] 

where Ki^{x,^p) and Piji^Xj^jp) are respectively the generalized photon-parton and purely 
partonic kernels. Usually, these equations can be diagonalized in terms of the singlet and 
non-singlet quark fragmentation functions as well as the gluon fragmentation function. In the 
following, we shall discuss the global features of the solutions of these evolution equations 
and perform an expansion in powers of as of these solutions. For this purpose, several 
simplifications can be consistently made. 

As we already mentioned previously, it turns out that the gluon-to-photon fragmentation 
function is by orders of magnitude smaller than the quark-to-photon fragmentation func- 
tions. As can be seen in Fig. ^ the inclusion of contributions involving the gluon-to-photon 
fragmentation function in the evaluation of the photon production cross sections leads only 
to negligible changes in the result, particularly at large x. It is therefore legitimate when 
first examining the evolution equations to ignore it and we will do so throughout this sec- 
tion. A further simplification is made by ignoring the possible transitions between different 
quark flavours, which yield a small contribution to alP^}}^. As a direct consequence of these 
considerations, the flavour singlet and non-singlet quark-to-photon fragmentation functions 
become equal to a unique fragmentation function Dg^^ which satisfies an evolution equation 
given by, 

dDg^^{x,iJ,j,) 



d\n{fi 



2 




(3.2) 



The generalized splitting functions Pj^ and Pij have a perturbative expansion in the strong 
coupling as{fi%). Explicit forms for the leading and next-to-leading splitting functions ap- 



pearing in these expansions can be found in |11]. In particular. 



Pg,{x,^il) = p!,%{x) + {^) p!,%{x) + ■ ■ ■ , 

PU^.f^l) = (^) + • • • ' (3-3) 
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with PW^^ pW^ and being the leading and next-to-leading order quark-to-photon and 
the leading order quark-to-quark splitting functions encountered in Section The resulting 
evolution eq. (|3.2|) has a similar form to the next-to-leading order evolution valid in the 



fixed order approach, eq. ( |2.9D . However, unlike in eq. (p.9| ), the strong coupling as is now 



a function of the factorization scale. As usual, the running of the strong coupling, a^, is 



determined by the beta function |jT4| 



where Po = (33 - 2Nf)/6 and pi = (306 - 38A^/)/6/?o- 

As in our fixed order approach, the full solution of the inhomogeneous evolution equation 
is given by the sum of two contributions; a pointlike (or perturbative) part -DgL^ which is 
a solution of the inhomogeneous equation ( p.2| ) and a hadronic (or non-perturbative) part 



D'^'^^ which is the solution of the corresponding homogeneous equation. 

Approximate solutions of the evolution equations are commonly obtained as follows P, 0] . 
First an analytic solution in moment space is obtained in the leading logarithm (LL) or 
beyond leading logarithm (BLL) approximations. These are then inverted numerically to 
give the fragmentation function in x-space. At LL only terms of the form (a" In"^^ are 
kept while at BLL both leading (a" In""*"^ ix^p) and subleading (a" In" logarithms of the 
mass factorization scale are resummed to all orders in the strong coupling a^. The strong 
coupling itself is obtained by integrating eq. ( p.4| ) and retaining only the first term in the 
LL case, while keeping both terms at BLL. We shall examine two approximate solutions of 
the evolution equation (|3.2|) more closely below. 



3.1.1 The pointlike part of D, 



Let us first concentrate on the pointlike part of the fragmentation function. The moments 
of the fragmentation function are defined as. 



•J 

In moment space, the leading logarithmic solution takes the form P, 



a[n] 



while the beyond leading logarithmic solution reads 0, p!5[| . 



27r 



27r 



as(/ii 



a{n) + h'{n) 



.«s(/io)/ 



(3.5) 



(3.6) 
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+b{n) 



(3.7) 



with Pj!^q(n) denoting the moments of the leading order quark-to-quark sphtting function. 

Independently of the precise definitions of the functions a, h and h' which we will come 
to next, both LL and BLL solutions have an asymptotic behaviour given by, 

,2^ 



27r 



2ti ) asifip) 



a{n). 



(3.8) 



This asymptotic form lends support to the common assumption that the quark-to-photon 
fragmentation function Dg^^ is O (a/as). This assumption is in contrast with that adopted 
in the fixed order approach (cf. Section 0) where the quark-to-photon fragmentation function 
is 0{a). It can lead to significant differences in the respective expressions of the one-photon 
production cross sections. We shall study these discrepancies more closely in Section ^ 

The functions a, h' and h are functions of the leading and next-to-leading quark-to-quark 
and quark-to-photon splitting functions in moment space. They are given by |p, |15[, 

1 P^%{n) 



a[n) 

b'{n) 
h{n) 



(n) 



/5o 



/3o 



/3o 



P!,\{n) - /?iP(^,(n)l a(n) 



P!,%in) 



{Pi'2 



{n 



n 



(3.9) 



Here Pa^f,(n) denote the moments of the mth order a ^ b splitting functions. They can be 
found in M- 



We will explicitly use these definitions in Sec. to estimate the difference between the 
numerical LL and BLL solutions in x-space obtained by inverting eqs. ( p.6| ) and ( p.7| ) p, 0] 
and analytic expressions obtained making a Taylor expansion in of these LL and BLL 
solutions. 



3.1.2 The hadronic part of D, 



9^7 



The hadronic part of the quark-to-photon fragmentation function is a solution of the homoge- 
neous evolution equation (eq. ( ^I2D with P^^ = 0). As for the solution of the inhomogeneous 
evolution equation, we can obtain solutions in the LL and BLL approximations defined 
above. In moment space we have. 



D 



had,LL , 
1^1 ' 




00 



D 



had,LL , 
q^-y I 



n,/io) + 0{as) 



(3.10) 
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and [lis 



D 



had,BLL , 
9^7 ' 



a 



D 



had,BLL , 



(3.11] 



To obtain the pointlike LL and BLL solutions of evolution equations (eqs. ( p.6|) and (|3.71 )), 
in any case, one needs to specify the non-perturbative input D^^^(n, /^o), which in both the 
conventional 0, |^ and fixed order approaches is proportional to a. Once this initial 
fragmentation function is chosen, within the conventional approach there are two different 
ways adopted in the literature to define the complete fragmentation function. 

One way is to consider the complete solution in a given approximation to be obtained as 
the sum of the pointlike and hadronic solutions within that approximation. In particular, the 
full BLL solution is obtained by adding the BLL pointlike and hadronic parts together. In 
this approach, which is adopted by Gliick, Reya and Vogt (GRV) in the full solution can 
be obtained by iteration. However, when solving the evolution equation in moment space, 
one inevitably includes terms which are beyond the order considered (as can be seen in 
eqs. ( p.lOp and (|3.11|) ). When the inversion into x-space is performed, one obtains spurious 
terms which can lead to significant contributions and have to be systematically omitted 0. 

An alternative way to construct the full BLL fragmentation function, adopted by Bourhis, 
Fontannaz and Guillet (BFG) in 0|, is to associate the BLL pointlike part with the LL 
hadronic part. From eq. ( p. 121) , it appears that the treatment of the hadronic part of the 
fragmentation function in this second conventional approach is conceptually closer to its 
treatment in the fixed order approach. 

At this stage it is important to mention that, in either of the two conventional approaches 
used to determine the quark-to-photon fragmentation function Dq^y{x, fip) described above, 
the hadronic input associated with the LL pointlike part turns out to be negligible and can 
be described by a VMD model P, . However, the hadronic input associated with the BLL 
pointlike part is sizeable and cannot be purely described by such a model anymore. We will 
come back to this important point and to the possible forms of this hadronic input in some 
detail in Section ^. 



3.2 Analytic expansion of ^^[^^{x, /ip) 

In this subsection we shall concentrate on the pointlike part of the fragmentation function 
and more precisely on obtaining an analytic expression for it by making a series expansion in 
as- As already mentioned, the LL and BLL resummed expressions in x-space of the pointlike 
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fragmentation function can be obtained by inverting numerically eqs. (|3.6| ) and (|3.7| ). Ap- 
proximations of these resummed solutions in x-space can however be obtained analytically. 
First, one expands the expressions for the resummed fragmentation functions in moment 
space as a series in as, up to a given order. The truncated series can then easily be in- 
verted analytically to yield an approximate expression for D^_^^ in x-space. These expanded 
expressions of the pointlike quark-to-photon fragmentation function can then be compared 
with the LL and BLL resummed expressions of -D^i^^ which are only known numerically. 

More precisely, an analytic expansion (up to 0{aas)) of the LL expression for the frag- 
mentation function is obtained as follows. First, the LL expression for as(/io) truncated at 
order reads. 



a 



{LL) 



1 - 



1 + 



2tt 




In' 




0{al 



(3.12) 



Inserting this expression into eq. 



and expanding in series up to order ag^iMp) one obtains. 





9^7 



n] 



(3.13) 



which can be trivially inverted to yield the expansion of the LL pointlike fragmentation 
function in x-space. 



+0(al). 



1 I ae 




(3.14) 



Similarly, an expansion (up to 0{aas)) of the BLL pointlike fragmentation function is 
obtained by considering the BLL expression of asifi^), i.e. the expression obtained retaining 
the term proportional to f3i, which reads. 
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(3.15) 



Inserting it into the resummed expression of the BLL fragmentation function given in 
eq. ( p.7|) , the expanded expression of the BLL pointlike fragmentation function in x-space 
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Figure 5: The pointlike quark-to-photon fragmentation function for the up-quark evaluated 
at /ii? = Mz at (a) LL and (b) BLL. In each case, we show the numerical resummed prediction 
as well as the difference with the corresponding fixed order expansions given in eqs. (|3.14|) 
and ( |3.16| ). The hadronic scale is (a) yUo = 0.50 GeV for LL GRV and (b) /iq = 0.55 GeV 
for BLL GRV and /xq = 0.71 GeV for BLL BFG. In (a) the difference between the LL GRV 
result and the leading term in the perturbative expansion eq. ( |3.17| ) is shown short-dashed. 



thus reads, 



^pZ,BLL(e.p.)(^^^^) - D^^^^^^^^^\x^^,) + (^) In (^^) P«^(x) + 0{al). (3.16) 



The expanded expressions for the pointlike part of the quark-to-photon fragmentation 
function can now be compared directly with the numerical solutions of the LL and BLL 
resummed expressions for the pointlike fragmentation function for a fixed value of fip and 
over the whole x range. These comparisons are shown in Fig. ^ using the parameterization 
of the LL and BLL fragmentation functions given by Gliick, Reya and Vogt (GRV) in 0] and 
the BLL parameterization given by Bourhis, Fontannaz and Guillet (BFG) in 0. In each 
case, we show the fragmentation function for the up-quark multiplied by the momentum 
fraction x. The other flavours have a similar behaviour. To compare the series expanded 
fragmentation function with the resummed expression, we also show the differences between 
the resummed and the expanded solutions given in eqs. ( |3.14|) and ( p.l6|) for the appropriate 
choices of fip and /iQ. That is /iq = 0.50 GeV for LL GRV, /iq = 0.55 GeV for BLL GRV and 
Hq = 0.71 GeV for BLL BFG. In all cases we choose hf = Mz- As a further comparison, in 
Fig. ^(a) we also show the difference between the resummed LL fragmentation function and 
the 0{a) term of the series expansion. 
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P{0) infe 
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(3.17) 
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Inspection of Fig. ||(a) suggests that D^j^^ which is the only term present in the lowest 
order solution obtained in the fixed order approach is insufficient to correctly reproduce the 
behaviour of the LL resummed expression. On the other hand, we see that in the region 
0.2 < a; < 0.9 the expanded expression of the fragmentation function is remarkably close 
to the resummed expression for both LL and BLL pointlike solutions and in the BLL case 
for the both parameterizations (GRV or BFG). For small x < 0.2, there are possible 
large logarithms of x in addition to contributions from the gluon fragmentation function 
that are not treated correctly in the expanded result. At large x, x > 0.9, there is also a 
significant difference between the two approaches. This discrepancy could be a sign that 
large resummation effects are present. Or it could indicate that the presently available 
parameterizations for the resummed fragmentation functions are not accurate at large x and 
particularly for x > 0.95. In fact, this discrepancy can be traced back to the presence of 
logarithms of (1— x) that are explicit in the expanded result. These logarithms should also be 
present in the numerical resummed results. However, the parameterizations are necessarily 
obtained by inverting only a finite number of moments and it is a well known problem to 
describe a logarithmic behaviour with a polynomial expansion. 

As the resummed fragmentation functions were obtained after the ln(/i|n) had been re- 
summed, the general agreement with the unsummed and expanded fragmentation func- 
tions leads us to question the necessity of such a resummation at LEP energies. Moreover, 
this agreement in the BLL case has another important consequence. As can be seen from 
eq. the expanded expression for the BLL pointlike quark-to-photon- fragmentation 



function is also equal to the expression of the next-to-leading perturbative fragmentation 
function obtained in the fixed order approach as given in eq. ( p.l2|) , where the hadronic 
input is neglected, D^^^{x, fio) = 0. It is therefore instructive to implement the expanded 
expression for the quark-to-photon fragmentation function in the evaluation of observable 
one photon cross sections at LEP energies. Indeed, doing so will enable us to compare the 
results for these cross sections obtained in different approaches and to isolate easier the 
differences between them, a task to which we will now turn in Section 0. 



4 The cross section in the different approaches com- 
pared 

We are finally interested in comparing the inclusive and 'photon' + 1 jet cross sections 
evaluated in the two essentially different approaches, at fixed order and following the con- 
ventional approach. In Section |^, we have described how the fragmentation function and 
the one-photon cross section are defined in the fixed order approach, while in Section ^ 
the derivation of the LL and BLL expressions of the fragmentation function has been dis- 
cussed. These resummed expressions were determined within the conventional approach, as 
approximations of the solution of an all-order evolution equation ( |3.2| ). By making a Taylor 
expansion in the strong coupling, analytic expressions for the LL and BLL pointlike solutions 
were also considered in Section |3.2| . Nothing however has been said so far concerning the ex- 
pressions for the single photon production cross section within this conventional formalism. 
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We shall fulfill this task in this section. 

In the following, we shall consider four different classes of expressions for the one-photon 
production cross section. These classes will be defined depending on whether the resummed 
(LL or BLL) expressions of the quark-to-photon fragmentation function or the expanded 
expression of the fragmentation function as given by eq. ( |3.16|) are used in the cross section. 



Secondly these classes will be determined depending on whether the direct contributions 
to the cross section are evaluated as a perturbative series in as up to 0{as), or whether 
these direct contributions are evaluated by using a conventional power counting, associating 
the powers of as and the powers of In together. The results obtained for the inclusive 
and 'photon' -|- 1 jet cross sections following any of these four approaches to evaluate the 
one-photon production cross section will be compared to the OPAL and ALEPH data in the 
forthcoming section. 

As we will see, in the category of approaches using the expanded expression of the 
quark-to-photon fragmentation function in the cross section, the specification of the input 
function D^l^^{x, fio) in a given factorization scheme is an important and subtle point which 
is treated differently by GRV in and by Bourhis, Fontannaz and Guillet BFG in [|]. We 
shall describe the determination of this input fragmentation function according to either of 
these two groups in some detail in the second part of this section. 



4.1 Approaches using the resummed Dq^j{x, fip) 
4.1.1 Direct contributions evaluated at fixed order in as 

Let us first concentrate on the expression of the one-photon production cross section in the 
MS factorization scheme obtained using the LL or BLL resummed expressions of the quark- 
to-photon fragmentation function while keeping the direct hard scattering terms building the 
cross section at fixed order in the strong coupling constant as up to 0{aas)- In this case the 
cross section takes the same form as in the fixed order approach described in Section 0, and 
is described by eq. (|2.3| ) (with replaced by x) at 0{a) and by eq. (^1^) at 0{aas). Rather 
than the fitted forms (eqs. ( p.l3| ) and ( |2.14| )), at 0{a) the LL fragmentation function should 
be considered in eq. ( |2.3|) , while at 0{aas) the BLL fragmentation function needs to be taken 
into account in eq. ( [2.4] ). Provided, the solution of the all order evolution equation can be 
accurately determined, the cross section evaluated following this approach is theoretically 
preferred as it is the most complete: It includes all direct terms up to order aas and all 
fragmentation contributions proportional to (a" In" fXp) and (a" In""*"^ Hp) at all orders. 

In order to evaluate either of the single photon production rates according to this prescrip- 
tion, we simply need to replace the fragmentation functions defined at fixed order appearing 
in these leading and next-to-leading cross sections by the resummed LL or BLL fragmenta- 
tion functions and leave the remaining terms in the cross sections unchanged. Consequently, 
by following this approach one is in principle able to test the universality of the LL and 
BLL fragmentation functions, in particular when these functions are employed to evaluate 
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the 'photon' + 1 jet rate, an observable which was not used to determine these functions. 



4.1.2 Direct contributions conventionally evaluated 

Second, let us consider the expression of the cross section in the MS scheme arising when 
one uses the resummed LL or BLL fragmentation functions but when one considers the 
direct terms of the single photon production cross section evaluated with the conventional 
power counting. These terms are obtained by keeping only the leading or beyond leading 
logarithmic terms of the factorization scale fip up to a given order in ct^. As it is commonly 
done 0, 1^, we shall follow this prescription, which defines the conventional approach to 
obtain the LL and BLL cross sections, and only retain terms proportional to (a" In""*"^ 
at LL and terms of the form (a^ In"'*'^ Hp) and (a" In" /i^) in the BLL expression of the cross 
section. 

Remembering furthermore that in this conventional approach, an inverse power of is 
associated with the quark-to-photon fragmentation function as in eq. (|3.8|) , the LL and BLL 
expressions of the cross section are simply, 

q=l 

E (D,^yi^,f^F) + (^) ® + (^^^ (x,/i^)^(4.1) 

As already mentioned before, due to the different power of associated with the fragmen- 
tation function in this approach and in the fixed order approach, the terms present in these 
two equations and in the corresponding fixed order equations (( p.3|) and (^^ )) differ sub- 
stantially. As explained at length in ||^, this conventional procedure of associating an inverse 
power of as with the fragmentation function is clearly appropriate when the logarithms of 
the factorization scale are the only potentially large logarithms but is theoretically less 
consistent when different classes of large logarithms can occur as in the 'photon' -|- 1 jet 
cross section. 



(To dx 

1 da^^^ 
(Jo da; 



4.2 Approaches using the expanded Dq^^{x^ fip) 

In this subsection, we shall consider the formulations of the cross section obtained using the 
expanded (up to 0{aas)) expressions of the pointlike and hadronic part of the fragmentation 
function Dq^y{x, fip). In particular, we will consider the pointlike part of the fragmenta- 
tion function given by the expanded expression of the resummed form for the BLL pointlike 
fragmentation function discussed in Section |3.2| and defined in eq. ( |3.16|) . An expanded 



expression of the hadronic part will be given below. Finally, we shall implement the ex- 
panded expression obtained for the sum of pointlike and hadronic parts of the fragmentation 
function in the NLO and BLL formulations of the cross section given by eqs.( |2.4| ) and ( ^.1| ) 
respectively. 
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4.2.1 Possible definitions of Dg^^{x, fio) 

In order to obtain a definite prediction for tfie cross section using eitlier tlie GRV or BFG 
models of the quark-to-photon fragmentation function, we need to know how the hadronic 
part of the fragmentation function behaves and in particular we need to know how the 
input function D^E^^{x, /iq) is defined. Since within this context we have an expanded (up to 
0{aas)) form for the pointlike part of the fragmentation function, we here choose to consider 
also an expanded expression (up to 0{aas)) for the hadronic part of this function, given by, 

^ja.,(e.p.)(^^ /i^) = D^"^^(x, /io) + (g) In (^^) p!,% ® D,Z,(a:, /io) + Oiaal). (4.2) 



This expression is obtained by expanding eq. ( |3.11D and takes exactly the same form as 



the non-perturbative part of the fragmentation function defined in the fixed order approach 
described in Section Recall that in the fixed order approach, the input fragmentation 
function D'^^^{x, fio) present as a boundary condition in eq. ( |4.2|) was determined at each 
order by comparing the fixed order 'photon' +1 jet cross section with the ALEPH data, with 
leading order and next-to-leading order expressions given in equations ( |2.13| ) and ( |2.14| ) 



respectively. Note also that the fitted input fragmentation function determined within this 
fixed order context is clearly non-negligible at any order in a^. 

In the approaches of GRV or BFG , the treatment of the input fragmentation function 
Dg^^{x, fio) is quite different. At LL both GRV and BFG agree that D^^^^x, fio) is negligible 
and can be described by a vector meson dominance model (VMD) as explained in and 
respectively. However at BLL and in the MS scheme, the input fragmentation function 
cannot be negligible due to the presence of the direct term C^°) (see eqs. ( p.3|) , (|2.4| ) or ( |4.1|) ) 
and cannot be described by a VMD input alone. Indeed, C^°) diverg 1 and would 

drive the cross section to unacceptable negative values if a VMD input alone is considered 
for the input fragmentation function. Note that the requirement that the cross section is 
positive led the authors in ^, || to consider a term proportional to P^^!^|^ln(l — x)"^ in the 
expression of D^E^^{x, iiq). To summarize the discussion, in any resummed or fixed order 
approach, as soon as the direct term C^°) enters the cross section, as it does in the MS 
factorization scheme, the input fragmentation function D^!^^{x, (Iq) must compensate the 
large x behaviour of C^°^ 

So far all the formulae given for the cross section were given in the MS factorization 
scheme. It is well known however, that the direct and fragmentation contributions are not 
unequivocally defined. For instance in a different factorization scheme, S a part of C^°^ 
(direct term) AC^°) \s can be absorbed in the fragmentation function yielding new functions 
Dq^.y{x, fip) \s and C^°^ \s {x,fiF), so that. 



D,^j{x, fip) \s= Dq^^ix, fip) + [j^j ACf \s (x) (4.3) 
\s (x,^^) = {x,f^F) - ACf Is (x). (4.4) 



and, 
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In fact only the combination of both direct and fragmentation contributions present in the 
physical cross section is factorization scheme invariant. The approaches adopted by GRV or 
BFG use this factorization scheme ambiguity to determine D^^^{x, ^q), albeit in a slightly 
different way. 

The GRV group choose to work within a factorization scheme (called DIS^) in which the 
fragmentation input D^^^{x, /iq) at both LL and BLL is simply given by a VMD contribution 
and is therefore negligible. Essentially, within this new scheme the troublesome part of the 
direct contribution C!^^^ appearing in eq. ([4.1|) is removed by absorbing it into the definition 
of the quark-to-photon fragmentation function Dq^^{x, fip). More precisely, in this scheme 
we have, 

AC^f \dis, (x) = P^%{x) In ((1 - x)x') - 2 (i-^) , (4.5) 

and, 

D-^_^.^{x, /io) \dis= D-^^{x, fio) Ims + ACf \nis, (x) ~ D^^^^^ . (4.6) 

This transformation holds for any quark flavour. Within this scheme, the direct term 
C^*^-* \dis^ is regular as a; 1. The evolution equation of the quark-to-photon fragmen- 
tation function in the S-scheme is modified as well, 

Pi% \s= Pi% Ims - (^) AC^f \s ®Pi%. (4.7) 

Invoking a perturbative stability argument, this group considers the evolution equation in 
the D/S'-y- scheme, solves this equation in that scheme and then since the MS-scheme is 
traditionally preferred in the evaluation of cross sections at higher orders, transforms back 
their results to the MS-scheme according to eq. ( [4.31 ). For illustration, the GRV up-quark 
fragmentation function a.t fip = Mz and in both the DIS^ and MS factorization schemes is 
shown in Fig. ^ 

Note that, to obtain the results in the DIS-y scheme, we have added a term proportional to 
\diSj (x) to the MS fragmentation function provided numerically by GRV. We see that 
the difference between the full and pointlike fragmentation functions due only to the VMD 
input is small in either scheme, and particularly at large x. Furthermore, we see that although 
the fragmentation function in one scheme may be well behaved as x — 1, in a different 
scheme, it will diverge as ln(l — x). Surprisingly, it appears well behaved in the MS-scheme 
rather than the specially constructed DIS^ scheme. As discussed earlier, this is due to an 
inaccuracy in the numerical resummed results produced by inverting only a finite number of 
moments. If the large x region is treated correctly, the DIS^ fragmentation function should 
be well behaved while the MS fragmentation function will exhibit a logarithmic enhancement. 

In order to be able to implement an expanded form for the complete quark-to-photon 
fragmentation function in the BLL and NLO (MS) expressions of the cross section given 
by eqs. ( |2.4| ), ( [4.1|) we need to know such a form in that scheme. The sum of the ex- 
panded expressions for the pointlike and hadronic parts of the fragmentation function given 
in eqs. ( |3.16| ) and ( [4.2| ) in the MS-scheme reads, 

D^^Z',^ Ims (a:, ^^F) = /^^^i^'^ Ims Z^^) + D'^tt""'-^ Ims f^r)- (4.8) 
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Figure 6: The GRV up quark-to-photon fragmentation function xDu^'y{x, fip) evaluated 
fip = Mz- The full (pointlike) fragmentation function is shown as a solid (dashed) line in 
the MS-scheme and as a short-dashed (dotted) line in the DlS^-scheme (see text). 



Rewriting the hadronic input in the DIS^ scheme yields, 

o2 



2tt J 

,2 



IDIS., [X) 




x,/io). (4.9) 



The last two lines of this equation are proportional to the VMD input, D^^^ \dis^, ix,fio), 
which is negligible 0], and can be neglected. Note that, an equivalent way to obtain this 
expression, with the model-dependent input -D^^^ Inis-y {x,fio) neglected, is to consider 
^5^^(a;,/io) Ims ^e given simply by -AC^ 1^^^^. 
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4.2.2 Direct contributions evaluated conventionally 

The expression for the BLL single-photon production cross section in the MS-scheme ob- 
tained using a conventional power counting of the direct terms and using the expanded 
expression of the fragmentation function given in eq. ([4.9|) (with D'^E^^ \dis-, (2;, /^o) = 0) 
then reads, 



4(17) 



+ (if) (S) (f ) « + (^) , (4.10) 

where we have only retained terms proportional to (a° In (aj In yu|.) and (aj In^ 



4.2.3 Direct contributions evaluated at fixed order in as 

The corresponding next-to-leading order expression for the cross section can however not 
be directly obtained by implementing the expanded quark-to-photon fragmentation function 



of eq. (^^ ) in the next-to-leading order MS cross section given in eq. (|2.4|) . We need to 
consider a further modification to the direct term present in eq. (p.4|) which is generated by 



the transformation of the non-perturbative input from the MS'-scheme to the DJS'.^-scheme, 
eq. (^4.61). With this change, the term proportional to Cf'> ®Dq^^{x, jjip) in eq. ( |4.1| ) becomes. 



+ 



'ae^\ /as 



27r / V27r/ 



"lew 




DIS^ 



Here, the term proportional to (7^°-* ® AC.^'') is genuinely of 0{as) and unlike the term 
proportional to Pjj'^^j could be ignored in the conventional BLL approach discussed above. 
However, for a consistent treatment at fixed NLO, it must be retained, so that, in the 
MS-scheme, 

1 d.-°(-. 1 d.«"<- . ^ /ae|X .a,. , ^ ^^^^^^ 



ctq dx (To dx 

Here, we clearly see which terms differ between the BLL and NLO cross sections when using 
an expanded expression for the fragmentation function. Another way to understand the 
origin of this additional term is gained by considering the expressions of the next-to-leading 
order cross section dcr^^*^ in both MS and DI S^-schemes. The requirement that these two 
quantities are equal, dictates that. 



'ael\ (as 



2tx / V27r/ ^ V27r 
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has to be invariant under the scheme transformation. Equivalently we have, 

C\'^ \dis,= C« Ims -C, ® ACf \nis, . (4.12) 

Let us now turn to the approach followed by BFG to determine Dq^^{x, fip)- Unlike 
the GRV group they consider the evolution of the quark-to-photon fragmentation function 
directly in the MS-scheme with a non-perturbative input given by, 

D^P_^^ Ims {x,f^o) ^ D^^'^^{x,f,o) - )(^), (4-13) 

where -D^^^ is fixed by a vector dominance model and is negligible. Effectively, the input 
is — C^™^\ the coUinear part of the direct term C!:^^ which is defined in eq. ( 2.16 ) and again 
diverges logarithmically clS X ^ -L. 

As mentioned before, the GRV group uses the same input and the same evolution equation 
for each quark flavour as we do in the fixed order approach. Unlike in our fixed order 
approach, where all 5 flavours are treated massless, the GRV group considers the masses of 
charm and bottom quarks {rric = 1.5 GeV and rrib = 4.5 GeV) and let the evolution equations 
of the heavy flavour fragmentation functions start at the appropriate mass thresholds 
In other words, for the charm-to-photon and the bottom-to-photon fragmentation function, 
they select Ho = 171^, rrif, respectively. The BFG group also considers the charm and bottom 
quarks to be massive. However, in their approach, the input distribution given in eq. (|4.13| ) 
is valid only for light quarks and the heavy quark input is treated slightly differently as 
explained in 0. For most applications though, (such as single photon production at LEP 
energies) we are far from the quark mass thresholds and the massless evolution of the charm 
and bottom fragmentation functions is a good approximation. For this reason we will not 
go into the details of the treatment of heavy quarks in the BFG approach and merely refer 
the reader to the original work 0. 

For reference, we show the full and pointlike up-quark fragmentation functions in Fig. ^ 
at /ii? = Mz- Here, the difference between the curves is not due to -D^^^ alone, but depends 
on the combination, D^^^^{x,no) — C^™'''(x). This should engender a significant difference 
at large x, but, because of the difficulty of obtaining accurate parameterizations at large x 
numerically, this has been obscured. 

The analogues of the BLL and NLO expanded expressions of the single photon produc- 
tion cross section following the BFG approach (at least for the light quark) are obtained by 
replacing AC!:^ \dis^ by in eqs. (|4.10| ) and ( [4.11j ). A similar transformation accompa- 
nied by replacing Dg^^{x,^o) |d/5^ by Z)^^^(x,/xo) in eq. ( |4.9D yields an expression for the 
MS quark-to-photon fragmentation function in the BFG approach. Finally, note that the 
only difference between the expression of the NLO single photon production cross section 
obtained following the fixed order approach with the BFG non-perturbative input given in 
eq. ( [4.13| ) and that obtained in the fixed order approach as described in Section ^ is the 
different non-perturbative input which in our approach is determined by the ALEPH data, 
see eq. ( p.l4|) . 



24 



0.02 



0.015 



full 

pointlike 



0.005 



0.2 



0.4 



0.6 



Figure 7: The BFG up quark-to-photon fragmentation function xDu^^/{x, fip) evaluated at 
fiF = Mz- The full (pointlike) fragmentation function is shown as a solid (dashed) line in 
the MS-scheme. 



4.3 Summary 

The different strategies for evaluating the single photon cross section described in this section 
can be summarized as follows; 

I) Resummed fragmentation function obtained numerically together with explicit power 
counting of the coupling constants, as described in sec. [4.1.1| and at NLO by eq. ( |2.4| )). 

II) Resummed fragmentation function with conventional power counting, i.e. associating 
an inverse power of as with Dg^^. See sec. [4.1.2| and eq. ( [4.1|) for the BLL result. 



Ill) Fragmentation function expanded as a series in as (given in eq. ([4.9|) and with the 
non-perturbative input D^t^^ Inis-y {x,^o) neglected) together with conventional power 
counting, i.e. associating an inverse power of with each power of In Hp. See sec. [1.2.2 
and eq. (g]TD|). 



IV) Expanded fragmentation function together with explicit power counting of the coupling 
constants, as described in sec. [4.2.3| and eq. ( [4.11|) . 



Provided the resummed solution of the all order evolution equation can be accurately 
determined, the approach using this solution and the direct terms evaluated at fixed order 
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(approach I) represents the theoretically preferred approach. The approach evaluating the 
direct terms at fixed order and using an expanded and thereby approximate expression of 
the fragmentation function has however important advantages. It enables an analytic de- 
termination of the fragmentation function and yields factorization scale independent results 
for the photon production cross section evaluated at a given order in as- Furthermore, the 
implementation of the non-perturbative quark-to-photon fragmentation function which was 
fitted to the ALEPH 'photon' + 1 jet in the fixed order expressions of the 'isolated' and 
inclusive rate yielded perturbative stable predictions which agreed well with the ALEPH 
and OPAL data. 



In the next section we shall see how the theoretical predictions obtained following any of 
these approaches and using the GRV or BFG schemes compare with the OPAL and ALEPH 
data. For definiteness, in approaches I and II we will use either the parameterization of the 
pointlike fragmentation function of the GRV group or the sum of pointlike and hadronic 
parts (set I) of the BFG group. In approaches HI and IV, we will consider eq. ( [4.9| ) with 
the VMD input -D^i^^ 1^/5^ (a;,/xo) set to zero to describe the evolution equations for all 
active fiavours in the GRV case, and the same equation with AC^°) Id/s^, replaced by C!^°^'^ 
in the BFG case. As a result, in these approaches and in either of the two schemes all 
quark fiavours satisfy a massless evolution equation. Finally, the light fiavours start their 
evolution at /io = 0.55 GeV and at /zq = 0.71 GeV respectively in the GRV or in the BFG 
schemes, while the heavy c and h quark fragmentation functions start to evolve at /zq = frtc 
and /io = respectively. 



5 Results 



In the previous sections we have completed a detailed comparison of the two fundamentally 
different approaches (fixed order and conventional) for computing photon cross sections. We 
have described how both the cross sections and fragmentation functions are defined in these 
two formalisms and in Section ^ we have given expressions for the cross section obtained 
using an expanded expression for the quark-to-photon fragmentation function Dq^^{x, fip). 
We have now collected all necessary ingredients to be able to evaluate the single photon 



production cross section in any of the four approaches summarized in Section [4.3| and could 
perform this task while defining the non-perturbative input according to either of the GRV 
or BFG groups. We shall apply these calculations to compute both the inclusive cross section 
and the 'photon' + 1 jet rate and make comparisons with the OPAL and ALEPH data in 
the second part of this section. 

Before turning to the cross sections, we first compare the analytic expanded expression 
of the quark-to-photon fragmentation function given by eq. ( [4.9| ) with the numerically re- 
summed BLL results for both GRV and BFG prescriptions. This is shown in Fig. |^ for the 
up-quark. Note that in addition to fixing the non- VMD input, AC^°) \s differently in the two 
schemes, the hadronic scale is also different, /iq = 0.55 GeV for BLL GRV and /iq = 0.71 GeV 
for BLL BFG. There are two ranges of interest, 0.2 < a; < 0.95 which is relevant for the 
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inclusive photon data from OPAL and 0.7 < x < 1 appropriate for the ALEPH 'photon' + 
1 jet data. 

We see that, except in the very high x region, the various fragmentation functions gen- 
erally agree well with each other in shape and magnitude. As discussed earlier, at large 
X, there are significant disagreements which are mainly due to deficiencies in the numeri- 
cal parameterizationsQ We therefore expect, that predictions for the inclusive photon cross 
sections (which run over a wide range of x) will be largely in agreement, while significant 
differences may be apparent in the 'photon' + 1 jet estimates which focus on the large x 
region. 

We also see that the expanded fragmentation functions defined according to the BFG and 
GRV prescriptions are quite different. This is in part due to the different choice of hadronic 
scale, but mainly due to the fact that the non-VMD BFG input is more negative than that 
for GRV. As can be seen from eq. ( [4.9| ) and the definitions of AC^^-' \Dis~f 

and C^™'), then 

for the same hadronic scale, 

^fT,)^^^ Ims p^f) - D^Z',^^^'' Ims ^^F) > 0. (5.1) 



5.1 The inclusive cross sections 

In this subsection we collect the results obtained evaluating the inclusive one-photon pro- 
duction cross sections following any of the four approaches described in Section 4. 

So far, we have ignored the gluon-to-photon fragmentation function throughout. To 
illustrate the tiny role the gluon-to-photon fragmentation plays in a physical cross section, 
Fig. ^ shows the BLL prediction for the inclusive rate for the pointlike GRV parameterization 
both with and without the gluon-to-photon fragmentation contribution. To make the small 
difference manifest, we have multiplied the gluon-to-photon fragmentation contribution by 
a factor of 100. We see that at large E^, even when multiplied by a factor of 100, the gluon 
fragmentation contribution is entirely negligible. At lower energies, the gluon fragmentation 
reduces the cross section by at most 5% at E^^ ~ 10 GeV. Similar results are also obtained at 
LL. In the following we shall therefore ignore the gluon-to-photon fragmentation contribution 
to the photon production cross section. 

5.1.1 Perturbative stability and /XiT'-dependence 

The results obtained using the GRV parameterization together with conventional power 
counting at LL and BLL (i.e. approach II) are shown in Fig. |T^. The upper/lower curves 
correspond to varying the factorization scale in the range 2M^/0.5M^. We see that the LL 

^In fact, it is striking to notice that the only quark-to-photon fragmentation functions which appear to 
diverge as a; ^ 1 are those which have at least one analytic component which itself diverges as a; — > 1. 
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Figure 8: The quark-to-photon fragmentation function xDu^^{x, fip) evaluated at fip = Mz 
in the (MS)-scheme. The NLO fit from the ALEPH 'photon' + 1 jet data is shown as 
sohd hne. The pointhke GRV (BFG) parameterization is shown dashed (dotted) while the 
expanded result of eq. is shown short-dashed (dot-dashed). 



and BLL predictions are similar and thus appear to be perturbatively stable. Furthermore, 
the factorization scale dependence is significantly reduced in going from LL to BLL. 

Let us see what happens when the inclusive rate is evaluated in approach I, i.e. using the 
resummed LL and BLL fragmentation function but with the direct contributions evaluated 



at fixed order. This is shown in Fig. |TT| at LO and NLO for the same three choices of 
the factorization scale as in Fig. [1^. In this case, it appears that we can draw the same 
conclusion regarding the fip dependence, it is reduced in going from LO to NLO. However, 
the LO and NLO results appear to be significantly different over the whole range of E^. This 
difference is caused by the presence of the direct term C..^^^ in the LO and NLO expressions 
of the cross section (see eq. ( ^.31 ) and (|2.4|)). The failure to describe the inclusive data with 
the LL resummed fragmentation function and with the direct contributions evaluated at 
lowest order in the expression of the cross section indicate that the direct and fragmentation 
contributions are not properly matched to each other when the cross section is evaluated in 
this way. 
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(GeV) 



Figure 9: The BLL prediction for the inclusive rate evaluated ai fip = Mz using the pointlike 
GRV parameterization with the gluon-to-photon fragmentation function ignored (solid) and 
multiplied by 100 (dashed). The experimental data is taken from p|. 



5.1.2 Comparison of results 



Fig. 12 shows the (NLO or BLL) inclusive cross section obtained using each of the four 
approaches described in Section ^ for iip = Mz- We show the predictions using both the 
GRV and BFG schemes while considering in each case the definitions of the fragmentation 
function given in Section 4. The various approaches give predictions which have a similar 
shape and lie in a common band which is well contained within the experimental error 
bars over the whole x range of the OPAL data. The agreement between the predictions is 
largely due to the similarity between the fragmentation functions as shown in Fig. ^ but 
also because the leading direct term C*^"-* is included in each approach. We note that the 
cross section obtained using the expanded expression for the fragmentation function (III and 
IV) lies by an almost constant amount above the prediction obtained using the resummed 
fragmentation functions (I and II) over the whole x range of interest. However, given the size 
of the experimental errors, all four predictions appear to describe the OPAL data equally 
well. 
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Figure 10: The (a) LL and (b) BLL predictions for the inclusive photon rate using the 
pointhke GRV parameterization and conventional power counting. The solid lines show the 
prediction for hf = Mz, while the short-dashed (long-dashed) lines show the expected rate 
for Hp = 2Mz (O.SM^). The experimental data is taken from 0. 



5.2 The 'photon' + 1 jet rates 



Let us now present the results obtained for the 'photon' + 1 jet rate. In the following, 
we focus on one particular value of the jet clustering parameter ?/cut, Vcut = 0.1. We note 
that in the range of interest for the ALEPH data, 0.7 < 2; < 1 the difference between the 
results obtained using the pointlike quark-to-photon fragmentation functions alone or the full 
fragmentation function (i.e. sum of pointlike and and hadronic parts) is small. Furthermore, 
the gluon-to-photon fragmentation plays an entirely negligible role. 



5.2.1 Perturbative stability and fip dependence 



Predictions for the 'photon' -|- 1 jet rate using the pointlike GRV parameterization together 



with conventional power counting at LL and BLL (i.e. approach II) are shown in Fig. 
We vary the factorization scale over a factor of 2 of the central scale fip = Mz, and, as 
before, the factorization scale dependence is significantly reduced in going from LL to BLL. 
However, the difference between the LL and BLL results is sizeable and the shape completely 
different. In particular, the BLL prediction does match the shape of the data quite well. 

Let us now analyze what happens if one considers the resummed expression for the GRV 
fragmentation function in an expression of the cross section where the direct terms are 



evaluated at fixed order in as (i.e. approach I). Fig. shows the LO and NLO 'photon' + 



1 jet predictions for the same three values of the factorization scale fip. As can be seen by 



30 




Figure 11: The (a) LO and (b) NLO predictions for the inclusive photon rate using the 
pointhke GRV parameterization and exphcit power counting. The sohd hues show the pre- 
diction for fip = Mz, while the short-dashed (long-dashed) lines show the expected rate for 
Hp = 2Mz (0.5M^). The experimental data is taken from 0. 



comparing the leading and next-to-leading order results, the factorization scale dependence 
is significantly reduced. We also see that the shapes of the curves displayed at leading and 
next-to-leading order are not dramatically changed. This could be viewed as an indication 
that the results obtained in this approach are perturbatively stable. For fip = Mz we find 
that, the lowest order prediction appears to be slightly below the data while the next-to- 
leading order result lies above the data. 



5.2.2 Comparison of results 



Fig. 15 shows the (NLO or BLL) 'photon' + 1 jet rates obtained using each of the four 
approaches described in Section H for up = Mz, with the fragmentation functions used in 
each approach as defined in Section 4. Ignoring the large z region where we have reason 
to doubt the accuracy of the parameterizations in methods I and II, we see that the BFG 
predictions lie systematically below that obtained using the GRV parameterization and go 
through the experimental data points. As discussed earher, this difference is due to both the 
choice of hadronic scale and the non-VMD input. The BFG input is smaller and the 'photon' 
-|- 1 jet data clearly selects this choice. Notice however, that the BFG parameterization for 
the fragmentation function unlike that of the GRV group was proposed well after the ALEPH 
data were released. As in the inclusive photon rate, predictions involving the expanded 
fragmentation function (approaches III and IV) always lie above the corresponding approach 
using the resummed fragmentation function (I and II). Again, the data clearly prefers the 
resummed fragmentation function. However, the shape of the predictions obtained with an 
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Figure 12: The inclusive photon spectrum estimated from the four approaches discussed in 
sec. ^ using (a) the GRV and (b) the BFG parameterizations. The experimental data is 
taken from 



expanded fragmentation function indicates that adding a negative constant to them would 
describe the data very well. 



6 Conclusions 



In this paper we have made a detailed study of photon production in hadronic events in 
electron-positron annihilation at LEP energies. First, we have used the fixed order approach 
of |lOl to estimate the inclusive photon spectrum and to compare it with the recent 



OPAL data Here, the fragmentation function is determined at large x {x > 0.7) by 
the ALEPH 'photon' + 1 jet data |l[] and is an exact solution of the evolution equation 
without resummation of logarithms of the factorization scale. As such, the prediction is 
scale independent and, surprisingly, agrees well with the OPAL data which corresponds to x 
values as small as 0.2. This is a powerful indication that the fragmentation function fitted 
to the ALEPH data is process independent and can be used to predict photon cross sections 
in other processes. 

Alternative methods to compute inclusive photon cross sections rely on numerically solv- 
ing the evolution equations with some non-perturbative input. This input has two pieces, a 
small vector meson dominance contribution together with a perturbative counterterm. Dif- 
ferent parameterizations deal with this ambiguity in different ways. We have examined the 



^Note that most of the events would be categorized as 'photon' + 1 jet events if a jet algorithm had 
been applied (See Fig. ||). 
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(a) LL 



- HF = Mz 
Hp = Mz 12 
tip = 2 Mz 



0.7 0.75 0.8 0.85 0.9 0.95 




Figure 13: The (a) LL and (b) BLL predictions for the 'photon' + 1 jet rate using the 
pointhke GRV parameterization and conventional power counting. The sohd hues show the 
prediction for hf = Mz, while the short-dashed (long-dashed) lines show the expected rate 
for Hp = 2Mz (O.SM^). The experimental data is taken from M. 



choices made by the BFG and GRV groups and used them to compute both the inclusive 
photon and 'photon' + 1 jet rates. To check the general behaviour of the fragmentation 
function, we have made an analytic series expansion in the strong coupling. As a result, 
we find that the large x behaviour of the fragmentation functions is not well reproduced by 
the parameterizations, the main problem being to describe a logarithmic behaviour with a 
polynomial. 

An additional subtlety is that although the fragmentation function appears to be 0{a), 
inspection of the evolution equation suggests a logarithmic growth with fip, and in many an- 
alyzes, it is ascribed a nominal power of a jag- Constructing the cross section at some partic- 
ular order depends on this assignment and different terms will contribute. This was discussed 
at length in Section ^. In addition, the gluon-to-photon fragmentation function is naively of 
0{aa^ and is expected to be much smaller than the quark contribution. This is indeed the 
case for physical cross sections in electron-positron annihilation where gluon production is 
suppressed, and we ignore the gluon fragmentation function contribution throughout. 

In order to better isolate the differences between the expressions of the cross section 
evaluated in a fixed order or in a conventional formalism, we have considered four ways of 
constructing the cross section for each parameterization; considering Dq^^ to be 0{a) or 
0(a/as) together with either the resummed or expanded solution of the evolution equa- 
tion. Provided the resummed solution of the all order evolution equation can be accurately 
determined, the approach using this solution and the direct terms evaluated at fixed or- 
der (approach I) represents the theoretically preferred approach. However, in Section ^ we 
pointed out that the approach evaluating the direct terms at fixed order and using an ex- 
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(a) LO 



- HF = Mz 
Hp = Mz 12 
tip = 2 Mz 



0.7 0.75 0.8 0.85 0.9 0.95 




Figure 14: The (a) LO and (b) NLO predictions for the 'photon' + 1 jet rate using the 
pointhke GRV parameterization and conventional power counting. The sohd hues show the 
prediction for hf = Mz, while the short-dashed (long-dashed) lines show the expected rate 
for Hp = 2Mz (O.SM^). The experimental data is taken from M. 



panded expression of the fragmentation function has some important advantages, such as 
eliminating the factorization scale dependence and having an analytic form. 

Predictions using these four approaches and either the GRV or BFG schemes were com- 
pared with the experimental data in Section |^. Reassuringly, in all cases, the NLO or BLL 
predictions are significantly less sensitive to the choice of factorization scale than the LO 
or LL predictions. We will therefore confine our comments to comparisons of the NLO or 
BLL predictions with the data. We note that estimates using the expanded fragmentation 
function systematically lie above those using the resummed fragmentation function. 

Unfortunately, photons can be confused with neutral pions. Consequently, the mea- 
sured photon cross section can only be obtained after a very large experimental background 
subtraction has been performed. As a result, the experimental errors are quite large. In 
particular, the OPAL inclusive photon data is unable to discriminate between any of the 
approaches or parameterizations used to predict the cross section at NLO or BLL. On the 
other hand, the ALEPH data does discriminate amongst the models, and, apart from the 
very high z region where the parameterization is suspect, prefers the BFG fragmentation 
function. Estimates based on either the expanded fragmentation function in both GRV or 
BFG schemes or the GRV parameterization for the resummed fragmentation function give re- 
sults that are systematically larger than the data allows. Nevertheless, the predictions based 
on the resummed BFG parameterization, with explicit power counting {Dq^^ of 0{a)) or 
conventional power counting {Dg^^ of 0{a/as)), agree well with the data but are too similar 
to be discriminated between. 
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(a) GRV 



(b) BFG 



0.7 0.75 0.8 0.85 0.9 0.95 1 



0.7 0.75 0.8 0.85 0.9 0.95 1 



Figure 15: The 'photon' + 1 jet rate estimated from the four approaches discussed in sec. ^ 
using (a) the GRV (b) the BFG inputs. The experimental data is taken from 



To summarize, we have shown how the inclusive and 'photon' + 1 jet data from LEP can 
be described by either the fragmentation function fitted to the ALEPH data or by the BFG 
solution of the evolution equation. In the latter case, the agreement needs however to be 
restricted to z-values below 0.95. We expect that this good agreement can be taken across 
to a variety of processes involving quarks and photons, such as prompt photon production 
at hadron colliders and the photon pair cross section at LHC. This may be of assistance in 
determining both the gluon content of the proton at moderate x values as well as in detecting 
a Standard Model Higgs-boson of intermediate mass via its two photon decay at the LHC. 
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